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ABSTRACT 

We present infrared interferometric imaging of the S-type Mira star^f Cygni. The object 
was observed at four different epochs in 2005-2006 with the IOTA optical interferometer (H 
band). Images show up to 40% variation in the stellar diameter, as well as significant changes 
in the limb darkening and stellar inhomogeneities. Model fitting gave precise time-dependent 
values of the stellar diameter, and reveals presence and displacement of a warm molecular 
layer. The star radius, corrected for limb darkening, has a mean value of 12.1 mas and shows a 
5.1 mas amplitude pulsation. Minimum diameter was observed at phase 0.94 ± 0.01. Maximum 
temperature was observed several days later at phase 1 .02 ± 0.02. We also show that combining 
the angular acceleration of the molecular layer with CO (Av = 3) radial velocity measurements 
yields a 5.9 ±1.5 mas parallax. The constant acceleration of the CO molecules - during 80% of 
the pulsation cycle - lead us to argument for a free-falling layer. The acceleration is compatible 
with a gravitational field produced by a 2.1+^ solar mass star. This last value is in agreement 
with fundamental mode pulsator models. We foresee increased development of techniques 
consisting in combining radial velocity with interferometric angular measurements, ultimately 
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allowing total mapping of the speed, density, and position of the diverse species in pulsation 
driven atmospheres. 

Subject headings: techniques: interferometric - techniques: radial velocities - stars: funda- 
mental parameters - stars: individual: x Cyg - star: AGB and post-AGB - infrared: stars 



1. Introduction 



Mira variables are low to intermediate mass AGB stars that pulsate with a period of about 1 year. They 
have a cool (T e ff < 3000^) and extended (R > 100/?©) photosphere. As suc h, they are br ight (M^ < -7) 



infrared beacons, individually observable far into galaxies of t he Local Group (IZijlstrall 19991) . They have the 
potential to probe places where the distance (eg, NGC 5128 in lRejkuball2004l) or reddening (eg, the Galactic 



Center in lMatsunaga et al.ll2009l) does not allow observation of the fainter/bluer - and rarer - Cepheids. 



However, the challenge to overcome is that Mira stars are both intrinsically complicated and ill- 
understood. Two important relations are of special interest: the period/luminosity (P/L) and the period/mass/radius 
(P/M/ R). The first relation has been derived from p opulation studies (sequenc e "C" in the LMC from lWood 
(12000 ) and also in the globular cluster 47 Tuc from Lebzelter & Wood (2005)). The present best parameter- 



ization of the P-L relation within our galaxy is (IWhitelock et al 



2008 



M k = -(3.51 ± 0.20)(log(P) - 2.38) - (7.25 ± 0.07) , 



(1) 



where P is the period in days. The zero point of this relation is the most uncertain parameter, with its 
dependence on the metallicity hardly known. The main difficulty is that parallax values are inaccurate and 
error-prone due to the large size and inhomogeneous surface brightness of the objects. 

The second relation, the P/M/R relation, has more relevance to the fundamental physics of the star. It is 
extremely dependent on the pulsation mode, but also, less crucially, on the surface density and metallicity. 
The P/M/R rel ation has bee n formally derived from numerical modeling of the fundamental pulsation mode 
of theses stars (|woodll 19891) : 



log(P) = -2.07 + 1.941og(7?//? ) - O.91og(M/M ), 



(2) 



Twenty years later, this model-derived relationship has still rarely been confronted with observation. This 
paper is a first step forward to establish the P/L and P/M/R relations on a new firm observational footing. 

Of the three crucial parameters (distance, mass, and radius), the angular diameter is far from being the 
easiest value to obtain. Because the surface gravity is several orders of magnitude lower than the sun, the 
pulsation of Mira variable leads to an extended atmosphere. In the cool upper layers, significant amount 
of the products of the helium fusion react to form di- and polyatomic molecules including TiO, SiO, CO 
and H2O. The forest of molecular lines and scattering from the dust lead to exotic intensity distributions 
not at all like a simple stellar disk. In the past, this substantially affected many stellar angular diameters 
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mode of pulsation (IBarfhes 



1998 



measurements (van Belle et al. 199 61: IPerrin et al.lll99 9). paving the way to contentious discussion on the 



Ya'Ari & Tuchman 



1999). 



However, now adays, interferometers are able to provide map s of the brightness distribution as a func- 
tion of wavelength (|Le Bouquin et al.ll2009l : IPluzhnik et al.ll2009r ). Images of the Mira star T Lep revealed 



a shell-like atmosphere, with a bright ch romatic zone distinctly detached from the photosphere. This could 
be the first image of what bhnakanboo^ called the MOLsphere, a zone of increased density in which for- 
mation of warm molec ular species would be favored. Accounting for this layer is th e key to obtain a correct 
value for the diameter (Perrin et al. I ll999l : Men nesson et ali 2002; Pe rrin et al. I l2004 . We will also show that 
we can apply to this layer a modified Baade-Wesselink method to derive the distance and mass of the star. 

The test star of this paper is x Cyg, a S-type Mira star. It has a pulsation period o f 408 days, a p hoto- 
metric magnitude ranging from 5.3 to 13.3, and intense emission lines at postmaximum (|Merrillll 19471) . This 
suggests a large pulsation amplitude. Images were obtained with the IOTA interferometer at four different 
stellar phases, chronologically cp = 0.93, 0.26, 0.69 and 0.79. In the next section, we describe the observa- 
tions and give a short overview of the dataset. In section [3j we use an image reconstruction algorithm to 
map the brightness distribution of the star. Precise geometrical parameters of the star, including existence of 
the molecular layer, are determined by model fitting in sectional From these values, temperatures and opac- 
ities are deduced in section [5] Finally, in section [6j we combine angular acceleration with radial velocities 
measurements to derive estimations of the distance and mass of the star. 



2. Observation and data reduction 



The interferometri c data presented h erein were obtained using the IOTA (Infrared-Optical Telescope 
Array) interferometer dTraub et al. 2003 ). a long baseline interferometer which operated at near-infrared 
wavelengths. It consisted of three 0.45-meter telescopes movable among 17 stations along two orthogonal 
linear arms. IOTA synthesized a total aperture size of 35 x 15 m, corresponding to an angular resolution of 
« 10 x 23 milliarcseconds at 1.65 fim. Squ ared visibility and closure phase measurements were obtained 
using the integrated optics combiner IONIC (|Berger et al.ll2003l ). IOTA ceased operation in July 2006. 



The declination of x Cyg made possible observations at low airmass with the IOTA interferometer 
during 6 months per year. With a period of 408 days, it allowed observations around both the minimum and 
the maximum brightness. It was observed at four different stellar phases: May 2005 {<p = 0.93), October 
2005 (<p = 0.26), April 2006 (0 = 0.69) and May 2006 (<p = 0.79). Observation information can be found 
in Table [H including dates and phase of observation, interferometer configurations and projected baseline 
lengths. Fig. Q] shows the u-v coverage achieved during these observation runs. The geometry of the IOTA 
interferometer and the position of the star on the sky constrained the extent of frequency coverage. 

The first three periods of observations were undertaken with a narrow H band filter (1.59//m< A < 
1.69//m). A dispersive setup was implemented shortly before our last ob servation run, resulti ng in a dataset 
featuring 7 spectral channels covering the range 1.5 /im< A < 1.8//m (IPedretti et al.l l2008). The science 
target observations are interleaved with identical observations of unresolved or partially resolved stars, used 
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to calibrate the interferometer's instru mental response. Th e in terleaved calibrator so urces (listed in Table [2]) 
were chosen in two different catalogs: iBorde et all (120021) and Me rand et al. (2006a). 



Reduction of the IONIC3 data was carried out using custom software (Monnier et al.ll2006r) . with ex- 
actly the same settings as the ones previously detailed in lLacour et al.1 (|2008l) . The output of this reduction 
softwar e are squared visib ilities (V ) and closure phases (CP). Data are available upon request in the OIFITS 
format (IPauls et alj|2005l) . They are presented in the 8 panels of Fig. [2] and [3] (superimposed on the data are 
curves of the best fit of the model presented in Sect. [4]). Eye analysis of the V 2 (upper panels) show a sig- 
nificant diameter variation between the different epochs. Also, at the longest baseline, the CP are clearly 
different from or n. This reveals the presence of an asymmetric brightness repartition, even though the 
complexity of CP estimators makes it hard to straightforwardly derive the level of asymmetry. 



3. Regularized imaging 



A first step in our data analysis was to convey the squared visibilities and closure phase into spatial 
information. The imaging was perfo rmed by the Mira reconstruction softw ar Jjl It stands for "Multi-aperture 
Image Reconstruction Algorithm" (IThiebautl |2008_|; iThiebaut et al.l 120031). For more details, we refer the 



reade rs to previous image reconstruction wor k using this software (ILacour et al 



2009) as well as more theoretical descriptions (ILe Besnerais et al.ll2008 



2008 



Le Bouquin et al 



Thiebaut & Giovannelli 



2009). 



Because the frequency coverage is far from being complete, the theoretical bijection between the fre- 
quency (Fig. [B and spatial domains is hard to achieve. It explains why the imaging algorithm requires a 
strong regularization term: the image is sought by minimizing a so-called cost function which is the sum of 
a regularization term plus data related terms. The data terms enforce the agreement of the model image with 
the different kind of measured data (squared visibilities and phase closures). The interpolation of missing 
frequency coverage is performed by the regularization and by strict constraints such as the positivity (which 
plays the role of a floating support constraint) and normalization. 

For x Cygni, the regularization term was a \ 2 minimization between the reconstructed image and a 
simple model of limb-darkened disk. Explicitly, the regularization term was: a £/t(4 _ nik) 2 where a is 
the relative weight of the regularization, is the intensity of k-th pixel in the sought image and is the 
flux of the k-th pixel given by the limb darkening model. The is a model of the brightness distribution 
of a stellar disk whose parameters (diameter and power of the limb-darkening) are beforehand adjusted on 
the squared visibilities. An advantage of this image reconstruction approach is that our prior favor radial 
symmetry. Hence this gives more strength to the relevancy of asymmetric features such as brighter regions 
in the restored images. 

The results are presented in Fig. 01 The four epochs are labeled by their stellar phase from upper left 
hand to lower right hand. The color scales of each image reconstruction is the same. The total brightness of 



1 http : //www- obs . univ- lyonl . f r/labo/perso/ eric . thiebaut/mira . html 
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each epoch is normalized according to the bolometric flux estimated in Sect. 15.1.31 

The variation in diameter is eye striking. Radius (defined by a 10% of the maximum surface flux) range 
from 9.4 (0 - 0.93) to 13.2 mas (<p = 0.26). This means a 40% expansion of the photosphere between the 
two epochs. Curves in Fig.[5]represent the radially averaged flux as a function of the radius. It is clear that, 
on top of diameter variations, changes in limb darkening are present, making it an important factor which 
has to be accounted for to determine precise photospheric radii. 



4. Model fitting 



4.1. Choosing the right geometrical description 



Interferometric observation of Mira stars are traditionally interpreted b y fitting simple brightness pro- 
files, such as a uniform disk (UD) or a Gaussian disk (GD) (lScholzll2003l) . Both fits showed inadequacy 
with our data. This is shown at the two first epochs of observation on Fig. [6j At low frequency, the data are 
best fitted with a Gaussian disk. At higher frequency, however, closure phases 7r-shifted imply the presence 
of a second lobe, ruling out a Gaussian type profile. With reduced^ 2 well over several hundreds, these fits 
confirm the need for a more complex geometrical description. 



Following an idea from lPerrin et al.1 (|1999T ). lMennesson et al.1 (|2002r ) first introduced a two-component 
model for Mira stars. It was a disk limb darkened by a close warm molecular layer. This model did well to 
account for spectral measurements in the K and L bands. However, we noted some problems when using 
this model to fit our measurements. First, because this model is spherically symmetric, it cannot account for 
closure phase measurements different from or 180 degrees. Secondly, the physical parameters (temper- 
ature and optical depth) present a degeneracy which cannot be resolved without spectral information. We 
therefore decided to use a purely geometrical description of the brightness distribution. This model contains 
seven parameters describing the presence of a photosphere with a center-to-limb darkening, a circumstellar 
envelope, and a spot on the stellar surface (Fig.|7]). 

The photosphere is modeled by a limb-d arkened disk. The center-to-limb variation (CLV) is a simple, 
single-parameter, power law: = fi a (|Michelson & Peasdll92ll : lHestrofferlll997|) : with ju = cos(#), 

and 6 the angle between the line of sight and the radial vector from the center of the star. Compared to a more 
classical 4-parameter law (as in Claret 2000), this law has the advantage of fitting various CLV shapes with 
a single parameter. The circumstellar envelope is modeled by an annular ring around the photo sphere. This 



is a si mple way to account for the existence of either (i) a warm molecula r layer as proposed b y lPerrin et al. 



(120041) or (ii) a molecular extension of the photosphere as modeled by llreland et all (|2008). Finally, we 
needed a way to account for the asymmetry. We used the simplest model available to explain an asymmetry: 
a single point-like spot. The model allows either a "cool" (dark) or "hot" (bright) spot. There is no restriction 
on the position of the spot (it is allowed to be outside the photosphere). 



To be sure that each one of these parameters are well constrained, it is crucial to know the influence 
they have on the visibility curve. The bottom panels of Fig.[6]give rough estimations of the zone of influence 
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of each parameter. At first order, the molecular layer is constrained by the low frequencies, the disk size by 
the first zero, the limb darkening by the height of the second lobe and the asymmetry by the closure phase. 
Of course, to a lesser extent, all the parameters affect each other at various degrees. 



4.2. Fitting the data 

An advantage of this model is that it consists in a sum of brightness distributions for which an analytical 
formula of the visibility function exists. The function of a power law limb-darkened disk, of parameter a 
and diameter 0*, writes: 

v, ^ y r(Qr/2 + 2) / -frvA) 2 \* ... 

V *M = g r(«/2 + * + 2)r ( * + i) [—^) ' (3) 

where v r is the radial spatial frequency (v r = Vm 2 + v 2 ) and T the Euler function. The visibility function of 
an annular ring of diameter #i ayer is: 

Vlayer(Vr) = Jo(27T9 layer V r ) , (4) 

where Jo is the Bessel function of the first kind. Finally, the visibility function of a hot spot writes: 

V spot (M, v) = exp {-liniX^u + F spo tv)) , (5) 

where the X spot and F spot are the coordinates of the spot (respectively right ascension and declination) relative 
to the center of symmetry of the stellar surface. The visibility of the full model is the weighted sum of the 
three visibility functions, hence: 

V(U, V) = F* V*(v r ) + Flayer Vlayer(Vr) + F spot V spot (w, v) . (6) 

F*, F^yer and F spot are the relative fluxes of, respectively, the disk, the envelope and the spot (F* + Fi ay e r + 

F sp ot = 1)- 

Closure phases are obtained by taking the argument of the product of three complex visibilities: arg {V{u\ , v\ ) • V(u2, V2) 
Squared visibilities are derived from the squared of the visibilities (\V(u,v)\ 2 ). Because u and v are wave- 
length dependent, accounting for the bandwidth smearing required to average the squared visibilities before 
fitting them to the data. For the 3 first observation runs, they are averaged over the 1.59- 1.69 /an bandpass. 
For the last run, the model was averaged over each channels bandpass (AA = 40 nm). The model is supposed 
to be achromatic. 

Thanks to the analytical expressions, the model can be swiftly fitted over a large range of parameter 
values. However, because the CP measurements are at the same time very precise (~ 1 degree) and very 
sparse (there is only one phase measurements for three squared visibilities), the spot position can have 
multiple^ 2 minima. To find the most likely set of parameters, we used the following strategy to achieve the 
global optimization of the x 2 '■ for a grid of given spot positions, we first map the x 2 minimized with respect 
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to the other parameters. Then we use the position which yields the best x 1 to initiate a local optimization 
with respect to all parameters by a Levenberg-Marquardt algorithm. 

The upper panels of Fig. [2]and[3]show the best fit of the model plotted on the squared visibilities (with a 
logarithmic scale). Since the asymmetry cannot be represented by a single radial visibility curve, we plotted 
three curves. The two dashed visibility curves are in the direction and at 90 degrees of the direction of the 
spot. The solid curve is the visibility curve toward the longest baseline measured. The difference between 
the model and the data points is plotted in the lower sub-panels. Residual errors on the V 2 average around 
1%. 

The lower panels of Fig. [2] and [3] show the best fit of the model plotted on the closure phases. Even 
though the point-like spot model is a very rough estimation of what the asymmetry could be, the general 
agreement between the model and the data confirms the validity of using such a simple representation. 
However, the results cannot exclude more complicated asymmetries, like the presence of multiple spots or 
other heterogeneities. For the October observations for example, small variations in the CP are not well 
matched by our model, hinting for a more complex repartition of the asymmetry. 

We noted that the minimum reduced^ 2 can be quite different from 1, ranging from a value of 1.3 
(March 2006) to 19.6 (May 2006). Lower x 2 can be obtained by using more complex models, like fitting 
two spots instead of one, or adjusting a chromatic limb darkening to the May 2006 data. However, the 
multiple parameters in such a case were too badly constrained to allow good determination. In the end, we 
decided to stay with our most simple model which gives a good compromise between fitting the data well 
and a reasonable number of free parameters. 

4.3. Parametric imaging 

The results and error bars of the fits can be found in Table [3] From these values, we can retrieve a 
brightness distribution of the object, ie, an image. The resulting images are presented in Fig. [U and should 
be compared to the regularized images of Fig. [4] Concordance is quite convincing, with a remarkable re- 
productivity in terms of diameter and position of the hot spots. The main difference is the presence of the 
molecular layer, situated around 0.5 stellar radius above the photosphere. The faintness of the layer (a faint- 
ness expected within our wavelength range of observation) may explain why the regularization algorithm 
did not image it. This on top of the fact that the regularization term does favor an empty environment. 
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5. Linking spatial to physical parameters 
5.1. Stellar diameter 

5.1.1. Rosseland radius 



In the case of an extended-atmosphere star such as Mira variables, a sensible definition of the radius 
as to be agreed upon. A quantity mostly used in Mira modeling is the Rosseland radius. This is an optical 
depth radius (Rr oss = r where tr oss = 1) and unfortunately not an observable quantity. In the case of the 
solar disk, the photospheric radius of the Sun is defined by the position of the CLV inflection point. 

The diameter values stated in Table [3] correspond to the furthest emission point of the photosphere, 
molecular layer notwithstanding. The problem of this definition of diameter is that it is highly model de- 
pendent. Such influence is emphasized in Table 0] The table compares the results of three different models 
fitted to the dataset of two epochs. It shows how profound the difference on diameter measurements can be 
depending on the model, well outside the range of error b ars. A solution co uld have been to consider only 
the uniform disk r adius, and to use tab les (for example in iDavis et alj|2000h to derive the Rosseland radius 
from simulations ([Ireland et al]|2004aUbl) . However, this technique poses the problem of fitting a uniform 
disk on a dataset which is no t compatible with it: the bias on the UD diameters depends on the Fourier 
coverage (IPaladini et al.ll2009h . 



However, a noticeable advantage of this work is the good coverage of spatial frequencies at all epochs. 
Our model, which covers a large range of CLV possibilities, allows to disentangle the main features of 
the stars (layer, spot). Simulations of pur e-continuum brightness profiles have steep flanks which mark 
the position of continuum-forming layers ( Scholz & Takedalll987l ). We are confident that the steep flank 
observed on the fitted model marks, as rightly as possible, the limit of the Rosseland Radius. 



5.1.2. Comparison with other interferometric observations 



Numerous x Cyg diameter measurements are present in the literature, but the multiple techniques of 
determination ma kes comparison di fficult. The first near-infrared interferometric observations of x Cyg 
were obtained by lYoung et all (120001) using the COAST instrument. In the 1.3 /mi continuum bandpass, 
they obtained a Gaussian FWHM of 13.9 + 0.8 mas at <p = 0.83. The 44% discrepancy with our diameter 
of 19.04 ± 0.09 mas at cp = 0.93 is difficult to exp lain by the difference of phase only. However it can 
be explained by the study of iHofmann et all (119981) . In their paper, they showed that the ra tio between 



contin uum radius and Gauss radius should be around 0.6, a value close to what is observed here. IPerrin et al. 



(120040 used the beam combiner FLUOR on IOTA to obtain narrow band measurements in the K-band. They 
fitted a more complete model including a uniform disk and a molecular layer with a wavelength dependent 
optical depth. They found an apparent stellar diameter of 21.10 ± 0.02 mas at (p = 0.24. This value is 
compatible with our measurements of 20.90 ± 0.12 mas obtained by fitting a model without limb darkening, 
but not with the diameter of 26.25 ± 0.08 mas measured when accounting for the CLV shape. So the 
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discrepancy is logically explained by the fact that they did not account for the limb-darkening. 



Another important diameter value still to be mentioned is the one obtained by lWeiner et all (120030 using 
the LSI heterodyne interferometer. They reported a diameter value of 39.38 ± 4.02 mas at <p = 0.51 in the 1 1 
fim continuum bandpass. Even though the phase of observation is different from our s , this value is far from 
being compatible. A possible interpretation could be the one proposed by Ohnaka (|2004k suggesting the 
presence of a warm, spectroscopically "hidden", H2O molecular layer. It is interesting to note that the 1 1 /mi 
stellar diameter m easurements match our measurements of the position of the molecular layer, something 
already noticed by IPerrin et al.1 (12004). 



Concerning the limb darkening, this paper offers the first measurements on a Mira star. But other 
observations exist on other types of stars. Merand et all (I2006bh report ed g ood fits with a = 16 toward 
Polaris and 6 Cepheid, two variable Cepheids. IWittkowski et all (|2006h and lLacour et al.1 ((2008) observed 
M giant stars and obtained respectively a = 0.24 ± 0.03 for Menkar and a = 0.258 ± 0.003 for Arcturus. 
In comparison, x Cygni's a values (greater than 1) look large. However, one has to consider the size and 
pulsation of the photosphe re of Mira va r iables . In that respect, our values are compatible with the CLV 
simulations presented in I Jacob & Scholzl (120021) . 



5.1.3. Effective temperature 



The effective temperature of a star depends upon its angular diameter and its bolometric flux, 



4 

"J ' ^Bol 



(7) 



where cr is the Stefan-Bol t zman n radiation constant, and Fb i is the observed flux, integrated over all wave- 
lengths. IWhitelock et all (|2000 ) reported J, H, K, and L band observations of the star at several phases. 
Thus, a phase-dependent bolometric flux can be estimated by fitting a black body distribution on the re- 
ported magnitudes. In Table [5] we present the photometric estimations as well as the bolometric flux at the 
epochs of observations. Applying Eq.|7j we deduced the effective temperatures (reported in Table [6b . They 
are interesting in the sense that they are much cooler t han what was deduced in IPerrin et aD (|2004l ) and are 
consistent with an M8 spectral type (IPerrin et al.lll998l) . 



5.1.4. A time-variable view of x Cygni 



Fig. [9]represents the flux, radius and temperature as a function of time. We fitted a sinusoidal model to 
the data. Best fit values are reported in Table [7J Linear radius is equal to 12.1 mas and the amplitude of the 
pulsation is 5.2 mas (43%). Minimum diameter happens at (p = 0.94 + 0.01. As expected, the temperature T 
is anti-correlated with the diameter, with a slight lag of 8% of the period (» 30 day s). This anti-correlation, 
as well as a similar lag, is also observed on Cepheids stars (lAndrievsky et al.ll2005l) . 



On the other hand, the bolometric flux is mostly correlated with the diameter, an indication that the 



-10- 



variation in the bolometric flux is dominated by the variation of the size of the object. The temperature plays 
a minor role on the bolometric flux, but making it phase shifted by 0. 14 in advance of phase compared to the 
radius. It is interesting to note that the visual magnitude, strongly affected by the molecular environment, is 
not in phase with the bolometric flux, but instead is in phase with the temperature. 



There are not many time-variable radius and temperature figures in the literature. lYoung et al.l ((2000) 
observed a variation with a maximum diameter at phase 0.6, but observations were done at a shorter wave- 
length (905 nm) which is affected by TiO absorption/emission. No noticeable diameter variation was ob- 
served in the less affected 1.3/rni bandpass. On the other hand, iThompson et al.1 (120021) observed a clear 
variation in diameter in the near infrared (K band) toward the Mira star R Tri Mira. They observed a 10% 
decrease in the UD diameter between phase 0.77 and 0.91. It agrees with the variation in radius observed 
here on^ Cyg, but unfortunately the time coverage on R Tri does not allow a more through-full comparison. 
Noteworthy is the work done by Woodruff et al. ( 20081 ): they derived time-dependant UD diameter for eight 
Mira variable, including^ Cygni. Minimum diameter is observed at phase « 0.7, a value inconsistent with 
our work. This could be explained by the fact that they fitted UD diameters which do not account for the 
molecular environment of the star. 



Least, a plot similar to the Fig. |9]is present in the Ph.D. Thesis of lStreckerl (119731) . He used bolometric 
flux and 3.5/mi observations of x Cyg to deduce a time- variable temperature and stellar diameter. In terms 
of variations, his results match our dataset well, with a minimum diameter around 0.9 and a minimum 
temperature at cf> = 0.5. However, compared to our work, he overestimated the size of the star by a factor 
1.5 (» 36 mas) and underestimated its effective temperature (» 2000 K). 



5.2. Temperature and opacity of the molecular layer 



The highly chromatic brightness of the layer observed by IPerrin et al.l (120041) immediately suggested 
the presence of molecules, mainly H2O and CO. Because of its molecular nature, it is extremely difficult to 
disentangle the relative effect of temperature and opacity on the brightness. 

To derive a first order estimat ion of the temperature, we used the assumption of a grey atmosphere. 
Hence, from lReid & Mentenl (I1997T) . the temperature writes: 



layer 



Ti(l- ^-(^ Mayer) 2 ) 



(8) 



From this first order estimation of the temperature, we can also deduce the optical depth of the molec- 
ular layer (ri ayer ) by using the flux conservation relation: 



^layer B(A, ^layer) ^layer 



1 - exp (-Tlayer) 



F+ B(^,r + ) 6> 2 exp (-Ti ayer ) 

where B(A, T) is the Planck function and A the wavelength. Hence: 



Tlayer = In 



1 + 



flayer B{A,T*)e\ 



Fi, B(A,T layei ) 6l yev j 



(9) 



(10) 



-11- 



The temperature and optical depth are reported in Table |6] With an effective temperature ranging from 
1750 K to 2 000 K, the layer is s ignificantly cooler than the excitation state of CO (Av = 3) molecules as 
observed by iHinkle et all (119821). On the other han d, the temperature agrees well with multi-wavelength 
interferometric observation from IPerrin et al.1 (120041) . The optical depth of the molecular layer is very low 
(t < 0.1), something we e xpected in the H band (and already reported on other Mira stars, as toward S Ori 
inlwitfkowski et al.1120081). 



6. Deriving Mass and Distance from the kinematics of the atmosphere 
6.1. On the existence of a molecular layer 



There is a strong debate about the p hysical nature of th e warm molecular layer as seen by optical inter- 
ferom etry. Precise interferometric values ( Perrin et alj2004 ) and im age reconstruction work ( Le Bouquin et al 
2009) tend to see a shell-like structure. Simulations (IScholal2003l) would favor a continuum emission from 
the photosphere up to a certain height in the atmosphere. This work does not pretend nor wish to resolve this 
issue. In the previous section, we chose to use the shell-like layer approach to understand our data, know- 
ing we are lacking the resolving power to convincingly distinguish between the two. From a geometrical 
point of view, the only difference between the two approaches only lies in the width of the layer, which is 
negligible in a shell-like representation. 

Whatever the width of t he molecular layer, clues on its k i nematics are offered by dynamical modeling 



of th e pulsating atmosphere (LB ertschin ger & Chevalier! 1 1985 



Bowen 



19881 : lBesselletal.il 19961) . Specifi 



cally, Bertschinger & Chevalierl (|l985l) showed that the atmosphere should be periodically criss-crossed by 
supersonic shocks. They computed a post-shock density 59 times the density o f the pre-shock, forming a 
zone in which rapid cooling would allow high nucleation rates ((Willsoru 120000 . The material is then in- 
wardly accelerated by the gravitational field, passing the sonic point closely below the shock-front. At this 
time, pressure forces become unimportant, making the rest of the trajectory ballistic-like. This theory nicely 
explains how warm molecules - and eventually grains - could be formed so low within the atmosphere of 
the star. 

If this model holds, the formation rates would decide the width of the molecular layer: density and tem- 
perature would define a spatially delimited zone in which molecules would happen to be more concentrated 
than in other atmospheric areas. 



6.2. A modified Baade-Wesselink method to derive the distance 



6.2.1. The CO (Av = 3) absorption features 



The 



Baadd (1 19261) and lWesselinkl dl946) method consisted originally in deriving the absolute diameter 



and distance of a star by means of photometry and spectroscopy. The interferometric Baade-Wesselink 
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method (applied to Cepheids in lKervella et al. b(X)4h differs by using a direct angular diameter measurement 
instead of photometry. We propose here to adapt this method to Mira stars, using the angular acceleration 
of the molecular shell. Concretely, the method allows deriving the distance by the relation: 



1 gangular (mas/s 2 ) 

parallax (mas) = — • 



P gradial velocity (AU/s ) 



(11) 



where p is a projection factor, ganguiar the geometric acceleration observed by interferometry, and g ra di a i velocity 
the absolute acceleration observed by radial velocity. An alternative method could have been to use the speed 
instead of the acceleration (first derivative of size instead of second derivative). The advantage of using the 
acceleration is to avoid a bias due to uncertainties in the local stellar velocity. 

Radial velocity of the molecular shell can be obtained through the CO second overton e vibration- 
rotation transitions (Av = 3), observable by spectroscopy around 1.6 fim (IHinkle et al.lll982h . The high 
vibrational energy of these molecules is characteristic of a warm environment (CO molecules in lower 
vibrational states are also observed but probe a cooler environment farther away from the photosphere). 
These excited molecules follow a roughly similar path on most of Mira stars, tightly correlated with the 
visual phase: they are created around maximu m visual brightness, steadily accelerated toward the star, and 



destroyed at the following stellar maximum (Hinklel Il978l : IHinkle et al.lll984l : iHinklel 1 1 995l : IHinkle et al. 



19971 : iLebzelter et al.ll2005f) . The excitation temperature is also correlated with the phase, showing a neat 
exponential cooling of the falling material. Line doubling is sometimes observed when the rotationally hot 
molecules form before complete dissociation of the les s excited ones. This behavior is making considerable 
sense in the light of a lBertschinger & Chevalier! (11985b scenario summarized in the previous section. 



One should be aware that using the Baade-Wesselink method in these conditions makes use of two 
strong hypotheses: 



The molecular layer expansion is radial. 

The molecular layer as seen by interferometry is the one from which originates the CO (Av 
absorption. 



3) 



The radial n ature of the CO disp lacement is a likely assumption because of the cycle to cycle repeatability as 
observed by IHinkle et all (1 19821) . The question of the concomitance of the molecular layer with the CO (Av = 



3) molecule is sensible, especially with respect to the unknown width of the molecular la yer. Regarding 
the effective temperature, the depth of the absorption lines reported by IHinkle et al.1 (119821) match 1700 K, 
close to the temperature of the layer reported in Table [6] On the other hand, the exc i tation temperature is 
higher, but the exponential decrease from 4400 K to 2200 K observed by IHinkle et al.1 (|l982T ) is hinting that 
excitation and effective temperatures are unlikely to be in equilibrium. 



In the near future, we expect high spectral resolution interferometers to lift this uncertainty by probing 
the apparent angular diameter of the star directly within the absorption line. 
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6.2.2. The projection factor 

The projection factor is one of the limiting parameters of the Baade-Wesselink methods. However, 
in the case of a Mira star imaged by interferometry, we can use the brightness distribution of the photo- 
sphere and combine it with the photosphere and molecular layer angular diameter. Assuming a shell of zero 
thickness, the projection factor writes: 

|f * 2nrp a dr 

- - (12) 



jf * Inr^ln^dr 

where ju* and ^i aye r are the projected radii, ie, = (1 - (r/R*) 2 ) 112 and yUi ayer = (1 - (r/Ri ayev ) 2 ) l/2 . Using 
this formula and the numbers stated in Table [6j the projection factor at the four phases of observations range 
from 1.22 to 1.27. In the following section we will use the projection factor p = 1.245 ± 0.025. 



6. 2. 3. x Cygni parallax 



In the upper panel of Fig.[l0]are plotted radial velocity data from lHinkle et all ([1982). For convenience, 
the data have been folded in phase, and shifted by 9.6 km/s. This velocity sh ift comes from an estimation 
of the local stellar velocity from CO measurements by mannier et al.1 (119901) . but is not relevant to derive 
the acceleration. We performed a linear fit over the period (0 < (p < 0.8) and derived a mean acceleration 
gradial velocity = -(1.10 ± 0.04) mrn.s 2 . Simultaneously, we fitted a parabola of constant inward acceleration 
over the position of the molecular layer (lower panel of Fig. [T0l . According to our measurements, the 
angular acceleration of the layer is ^angular — —(5.4 + 1.4) X 10 ^ mas.s 2 . Using Ecj. (|1 lb . this gives a 
parallax of 5.9 ±1.5 mas. The main source of uncertainty is coming from the estimation of the geometric 
acceleration of the layer, ganguiar- 

Several values of x Cyg's parallax exist in the literature, but no consensus exists yet, due to the difficulty 
of measuring the position of such a star of high angular diameter and inhomogeneous surfac e brightness . 
We can cite, among others, a measurement by the Allegheny Observatory (8.8 ± 1.9 mas, in 
and m ultiple parallax calcul ations obtained from the Hipparcos dataset: 9.43 + 1 .36 ma s in 



Stein 



1991) 



Perrvman et al. 



(| 19971) . 6.71 ± 1.00 mas in iKnapp et all (|2003l) . and 5.71 ± 1.12 mas in IVan Leeuwenl (|2007l ). Our mea- 
surement is compatible with this most recent value. Another way to measure the distance is to use the 
Period/Luminosity relation of Mira variables. According to Eq. (Q~|), the absolute magnitude ofx Cyg should 
be M k = -8.06 + 0.09. A ssuming a K magnitude of - 1 .95 (including interstellar or circumstellar extinction 
from IKnapp et al.ll2003l) . the distance of x Cyg should be 6.1 mas. This last value is also compatible with 
our measurement. 
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6.3. Mass derivation 

6.3.1. The gravitational field 
Assuming a free-falling molecular layer, the mass of the star can be derived by: 

M* = (13) 

where g = g ra diai velocity / P, G the universal gravitational constant, and /?i ay er the position of the molecular 
layer. Over the period < (p < 0.8, the average radius of the molecular layer is 18 mas, which gives 
^layer = 3-0+*g AU (for a parallax of 5.9 ±1.5 mas). Assuming a gravitational acceleration of g - -(1.37 ± 
0.05) mm.s -2 , this leads to a stellar mass = 2.1^^ M Q . Note that the uncertainty on the parallax 
dominates the error bar. 

We are fully aware that the free-fall approximation is an important assumption. Three other forces can 
affect the trajectory of the molecular layer: 

• pressure force 

• radiative pressure 

• centripetal force due to the rotation of the star 

The question is whether we can neglect them compared to the gravitational force. If not, our estimate would 
only yield a lower limit of the mass. 

Centripetal force could be dismissed due to the long rotation period of this evolved star. From conser- 
vation of momentum energy, it is likely that the star has a period well over several tens of years. 

Radiative pressure is responsible for the large mass loss. However, because of the small cross-section 
between CO and starlight, the kinetic energy transferred to the molecules would not result in a sensible 
acceleration. In other words, we expect /?co(Av=3) tne rat i° of radiation pressure to stellar gravity to be 
significantly below 1 . A dragging effect could be envisioned, but the velocity gradient between grains and 
molecules is necessarily small so low in the atmosphere. 



Gas pressure force in-between two shocks is likely to be small (Willson 2000). This can be under- 
stood as follo ws: the density gradient of simulated static atmosphere of Mira stars are extremely steep 
n a static atmosphere, pressure force induced by density gradient are rigorously equal to the 



gravitational force. In pulsating Mira stars, the atmosphere is "puffed up", with a density gradient several 
orders of magnitude smaller. Likewise, pressure forces will be several orders of magnitude lower than the 
gravitational force (up to a certain distance from the star). 

From an observational point of view, this idea is mostly confirmed by the constant acceleration of 
the CO (Av = 3) molecules. Pressure forces and radiative pressure cannot be constant over the pulsation 
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cycle. Bolometric flux varies by at least a factor 2, and pressur e forces may var y even more especially in the 
case of multiple ballistic trajectories per cycle as proposed by IWillsonl (120061) . Thus, if theses forces were 
dominant, we wo uld detect a non-lin ear effect on the acceleration of the molecules. Since the radial velocity 
measurements of Hinkle et al. ( 1982b rule out a difference in acceleration of more than 10% between phase 
and phase 0.8, it is unlikely that pressures forces slow the inward acceleration by more than a few percents. 



6.3.2. The Period/Mass/Radius relation 



Indirect confirmation of x Cyg mass is difficult. Indeed, 2.1 M© is quite massive for a Mira star 
(IVassiliadis & Woodlll993T) . However, x Cyg may not be a Mira of the m ost common sort. First, the star 
has often been noticed by the remarkable strength of his emission lines (eg Merrill 1947 ). Secondly , its 408 
day period is slightly longer then the average pulsation period of Mira stars (|Whitelock et al.ll2000h . Lastly, 



according to our parallax measurement, its mean radius of 440^° (12.1 mas) makes x Cygni a quite 
huge object. 



The feeling that x Cygni can be a quite massive Mira star is confirmed by the P/M/R relation. Indeed, 



applying x Cyg's parameters to Eq. (0 gives a mass of 3.1 
measurement. 



+2.7 
-1.2 



Mq. This value is compatible with our 



6.4. Discussion 

This work makes use of several assumptions which will certainly need to be fine-tuned. If anything, 
they emphasize the need to complete the method with detailed numerical modeling to assess the influence of 
atmospheric pressure forces. In the case of x Cyg, the consistency of the acceleration was the main argument 
to neglect pressure effects. We also neglected the influence of an eventual centrifugal force as it is unlikely 
because of the slow rotation of the star's atmosphere. Finally, we assumed the molecular laye r detected in 



the in terferometric data to be at the same position as the CO ( Av = 3) molecules observed by lllinkle et al. 



dl982h . This last point may be the most arguable one. Even though both spectroscopic and interferometric 
observations have been made in the H band, it is true that contamination may exist due to molecules (e.g. 
H2O) possibly formed at a slightly different height inside the atmosphere. Clear distinction between several 
possible molecular layers will nevertheless be possible in the near future thanks to high resolution spectro- 
interferometry. By probing the apparent stellar diameter within the depth of saturated absorption lines, it 
would provide the velocity and the position for individual species inside the atmosphere. This is a prospect 
which would not only strengthen our mass and distance determination method, but also would be bound to 
revolutionize our understanding of shock driven atmospheres. 

Evolved stars could prove to be precious secondary distance indicators. Compared to Cepheid, the 
fact that these stars are large and extremely bright in the infrared (where interstellar absorption is minimal), 
would make them a precious secondary indicator. To allow that, the P/L relation of LPV stars should 
be properly assessed. This work is a first step in that direction. Moreover, this method has an intrinsic 
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advantage compared to the traditional Baade-Wesselink method since the projection factor can be derived 
from the position of the molecular layer (even though we agree that the projection factor is not yet the 
limiting term in the final accuracy of this method). 

Deriving the mass of these stars also opens an important prospect to study stellar structure, chemistry, 
and mass loss. It is one of the main parameters (with radius) to constrain stellar instabilities and oscillation 
modes. This novel way to derive the mass is applicable to a wide variety of pulsating objects. The require- 
ment is to be able to probe a portion of the atmosphere where the gravitational field is significantly stronger 
than the pressure forces (conventional and radiative). It only happens in shock driven atmospheres, but nev- 
ertheless should be observable in numerous types of pulsating stars. For example, rece nt results on Cepheid s 
showed a faint upper atmospheric molecular layer likely to be gravitationally unstable (Merand et al. Il2006bh . 



Note that this technique does not require a regular pulsating variable; it can also be applied to any irregular 
or semiregular variables, featuring shock driven molecular layers. 

We acknowledge with thanks the variable star observations from the AAVSO International Database 
contributed by observers worldwide and used in this research. 
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Fig. 1 . — u-v coverage obtained on x Cyg at the four epochs of observation. Coordinates are in factors of 
the wavelength, ie, 10M/1 correspond to a baseline of 16 meters at A = 1.6 /an. Baseline lengths range from 
5 meters to 38 meters. 
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Fig. 2. — Squared visibilities (upper panels) and closure phases (lower panels) at the two first epochs of 
observation. The curves correspond to the best fit of the model described in Fig. [7] Three curves are plotted 
on the visibilities: the dashed curves are in the direction of the spot, and at 90 degrees from it. The solid 
curve is the visibility curve toward the longest baseline measured. The lower sub-panels show the residual 
errors. 
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Fig. 3. — Squared visibilities (upper panels) and closure phases (lower panels) at the two last epochs of 
observation. The curves correspond to the best fit of the model described in Fig. [7] Three curves are present: 
the dashed visibility curves are in the direction of the spot, and at 90 degrees from it. The solid curve is the 
visibility curve toward the longest baseline measured. The lower sub-panels show the residual errors. 
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Fig. 4. — Regularized imaging of x Cyg. The four epochs are labeled by their stellar phase, from upper 
left hand to lower right hand. The variation in diameter is eye striking with diameter changes of up to 40% 
between phase 0.93 and 0.26. Changes in limb darkening are also present. Hot cells on the photosphere are 
apparent, interestingly with a higher contrast of the cells when the star is at smaller diameter. The relative 
brightnesses are normalized according to H band magnitude from Table [5] The angular resolution of the 
IOTA interferometer is 10 X 23 mas (1.65/im). At% Cygni's distance of 170pc, 5.9 mas correspond to 1 AU 
(Sect. l6T3l . 
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Fig. 5. — Center-to-limb variation (CLV) at the four epochs of Fig|4j The flux is radially averaged. It shows 
a clear variation in the brightness distribution as a function of time. 
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Fig. 6. — Data of May 2005 (upper panels) and October 2005 (lower panels). On top of the measurements 
are plotted the best fit of a uniform disk and a Gaussian disk. While the UD has trouble to fit the low 
frequency measurements, the GD is unable to account for the CP 7r-shift. In the two lower panels are also 
summarized the "zone of influence" of the parameters used in the model presented in Fig. |7J the molecular 
layer (of flux Fj and radius /?/) will impact the data at low frequency, the diameter of the star (/?*) will 
determine the position of the first null, and the limb darkening (ji a ) will impact the higher frequency by 
modifying the height of the second lob. The asymmetry is clearly revealed by the CP. 



Fig. 7. — The model used in Sect. 0] It consists of a limb-darkened disk, a spot and a molecular layer 
represented by a ring. The free parameters are: 6+, #i ay er> en, X s , Y s , plus the brightness of the layer and the 
spot relative to the total flux. 
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Fig. 8. — Parametric imaging. The parametric model is showed in Fig. [7] The values used are the ones re- 
ported in Table [3] Apparency of the spot was chosen to be a fifth of the stellar diameter, below the resolution 
power of the interferometer. The relative brightnesses are normalized according to H band magnitude from 
Table|5J At^ Cygni's distance of 170 pc, 5.9 mas correspond to 1 AU (Sect. l(x2~3T ). 
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Fig. 9. — Temporal evolution of bolometric flux, radius and temperature of the photosphere of x Cyg- 
Values (Table © are plotted twice with +1 phase shift (white and black dots). A sinusoidal is fitted to the 
data points, giving a linear radius of 12. 1 mas and an average effective temperature of 2560 K (Table|7]). The 
temperature was derived according to: T e g = (4crF^ i/6^) 1 ^. It follows a periodic variation phase-shifted 
by 0.58 compared to the radius, ie, almost opposed to the dilatation of the photosphere. 
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Fig. 10. — Upper panel: Radial v elocity of the CO (A v = 3) molecule from iHinkle et all (119821) . folded 
in phase and shifted by 9.6 km/s (jwannier et al. 1990 ). The tilted lines result from a linear fit over the 
< 4> < 0.8 period. It corresponds to a constant acceleration of g rat ji a i velocity = -(1.10 ± 0.04) mm.s . 
Lower panel: ballistic trajectory of inward acceleration fitted to the position of the molecular layer. The 
highest point of the ballistic trajectory is at = 0.36. 
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Table 1 . Observation Log 



Date (UT) 





Configuration a 


Length b [m] 


2005 May 25 


0.92 


A15-B15-C0 


21 


2005 May 27 


0.93 


A15-B15-C5 


21 


2005 May 31 


0.93 


A25-B15-C10 


29 


2005 Jun 1 


0.94 


A35-B15-C10 


37 


2005 Oct 5, 6, 7 


0.25 


A5-B5-C0 


7 


2005 Oct 8, 9 


0.26 


A5-B15-C0 


14 


2005 Oct 10, 11 


0.26 


A15-B15-C0 


21 


2005 Oct 12, 13 


0.27 


A25-B15-C0 


29 


2006 Mar 29,31 


0.68 


A15-B15-C0 


21 


2006 Apr 2 


0.69 


A5-B5-C0 


7 


2006 Apr 7 


0.70 


A30-B15-C0 


33 


2006 May 11 


0.78 


A15-B5-C10 


15 


2006 May 12 


0.78 


A15-B5-C0 


15 


2006 May 13 


0.79 


A15-B15-C0 


21 


2006 May 14 


0.79 


A30-B15-C0 


33 


2006 May 15 


0.79 


A35-B15-C21 


38 


2006 May 16 


0.79 


A35-B15-C25 


38 



a Interferometer configuration refers to the location, in 
meters, of telescopes A, B, and C on the North/East, 
South/East and North/East arms, respectively 

b Length of maximum projected baseline 
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Table 2. Calibrators 



Calibrator Spectral Type UD diameter 



HD 176670 


K2.5 III 


2.330 


+ 


0.026 


HD 180450 


MO III 


2.770 


+ 


0.032 


HD 186619 


MOIIIab 


2.190 


± 


0.025 


HD 188149 


K4III 


1.490 


± 


0.020 


HD 197989 


KOIII 


4.440 


± 


0.048 
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Table 3. Best-fit parameters 





May/June 2005 


October 2005 


March/April 2006 


May 2006 




(0 - 0.93) 


(<p = 0.26) 


(<f> - 0.69) 


(<p = 0.79) 


0* (mas) 


19.04 ± 0.09 


26.25 ± 0.08 


23.97 ± 0.80 


21.49 ±0.11 


LD [a] 


1.34 ±0.05 


1.08 ±0.03 


2.54 ± 0.39 


2.42 ± 0.05 


flayer (mas) 


31.83 ±0.15 


40.75 ± 0.37 


35.48 ± 0.40 


27.35 ±0.13 


^layer/^total (%) 


6.5 ±0.2 


4.7 ± 0.2 


8.8 ±0.3 


8.13 ±0.2 


X spot (mas) 


5.22 ± 0.05 


-8.92 ± 0.39 


2.22 ± 0.42 


3.21 ± 0.20 


y spot (mas) 


2.97 ± 0.05 


2.96 ±0.10 


-4.24 ± 0.34 


-6.70 ± 0.09 


Fspot/F tota i (%) 


5.9 ±0.1 


1.7 ±0.1 


3.7 ±0.3 


1.7 ±0.1 


Reduced ^ 2 


6.8 


10.3 


1.33 


19.6 
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Table 4. Model dependent diameter measurements 







May 2005 (<p = 0.93) 




October 2005 (<p = 


0.26) 


Model 


UD 


UD + layer + spot 


LD + layer + spot 


UD 


UD + layer + spot 


LD + layer + spot 


(mas) 


16.24 ±0.07 


14.99 ± 0.07 


19.04 ± 0.09 


22.99 ±0.11 


20.90 ±0.12 


26.25 ± 0.08 


a 






1.35 ±0.05 






1.08 ±0.03 


flayer (mas) 




21.76 ±0.27 


31.85 ±0.15 




26.76 ± 0.24 


40.75 ± 0.37 


-flayer /-fTOTAL (%) 




10.9 ± 0.3 


6.5 ± 0.2 




15.1 ±0.6 


4.7 ±0.2 


X spo , (mas) 




5.47 ±0.14 


5.22 ±0.05 




-8.4 ±0.78 


-8.92 ±0.39 


r spot (mas) 




3.26 ±0.11 


2.97 ± 0.05 




2.15 ±0.17 


2.96 ±0.10 


F spot IF TOTAL (%) 




5.9 ± 0.2 


5.9 ±0.1 




1.2 ±0.1 


1.7 ±0.1 


Reduced x 


550 


29 


6 


536 


14 


10 
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Table 5. Bolometric Flux 



4> 


0.96 


0.26 


0.69 


0.79 


JD 


2453518 


2453653 


2453826 


2453869 


J (mag) a 


0.00 ±0.15 


-0.46 ±0.15 


0.15 ±0.15 


0.07 ±0.15 


H (mag) a 


-1.00 ±0.15 


-1.65 ±0.15 


-1.05 ±0.15 


-1.01 ±0.15 


K (mag) a 


-1.65 ±0.15 


-2.24 ±0.15 


-1.73 ±0.15 


-1.65 ±0.15 


L (mag) a 


-2.50 ±0.15 


-2.84 ±0.15 


-2.51 ±0.15 


-2.48 ±0.15 


^Bol 










(10" 13 Wcm- 2 ) 


6.83 ±0.38 


10.15 ±0.57 


6.80 ±0.35 


6.74 ± 0.36 



a Magnitudes obtained from iWhitelock et all (120001) 
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Table 6. Time dependent parameters of x Cyg 





May/June 2005 


October 2005 


March/April 2006 


May 2006 




0.93 


0.26 


0.69 


0.79 


0* (mas) 


19.04 ± 0.09 


26.25 ± 0.08 


23.97 ± 0.80 


21.49 ±0.11 


R*/Ro a 


348 ± 94 


480 ± 130 


439 ± 119 


393 ± 106 


r*(K) 


2742 ± 45 


2578 ± 40 


2441 ± 72 


2572 ±41 


flayer (mas) 


31.83 ±0.15 


40.75 ± 0.37 


35.48 ± 0.40 


27.35 ±0.13 


^layer /R-k 


1.67 ±0.02 


1.55 ±0.02 


1.48 ±0.07 


1.27 ±0.01 


^layer (K) 


1821 ±29 


1795 ± 28 


1747 ± 52 


2032 ± 32 


^"layer 


0.042 ± 0.002 


0.032 ± 0.002 


0.067 ± 0.007 


0.074 ± 0.002 




1.22 ±0.01 


1.26 ±0.01 


1.22 ± 0.02 


1.27 ±0.01 



a Assuming a parallax of 5.9 ±1.5 mas (Sect. l6T23T ) 
b Projection factor (Sect. 16.2.21 ) 
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Table 7. Sinusoidal fit a 



a b 0o 

F Bo i (10~ 13 Wcnr 2 ) 8.3 ±0.3 1.8 ±0.4 0.30 ± 0.03 

Diameter (mas) 24.2 ±0.1 5.1 ±0.1 0.44 ± 0.01 

Temperature (K) 2560 ± 30 226 ± 56 0.02 ± 0.02 



Sinusoidal defined by: a + b sin(27r(0 - 0o)) 



